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Remark on the perturbative component of inclusive r-decay 
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In the context of the inclusive r-decay, we analyze various forms of perturbative expansions which 
have appeared as modifications of the original perturbative series. We argue that analytic pertur- 
bation theory, which combines renormalization-group invariance and Q^-analyticity, has significant 
merits favoring its use to describe the perturbative component of r-decay. 
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A perturbative approximation in quantum chromody- 
namics as a rule cannot be exhaustive in the low energy 
region of a few GeV and a nonperturbative component 
has to be included. The reliability of extracting non- 
perturbative parameters from data is connected with un- 
certainties in the perturbative description of a process 
arising from the inevitable truncation of the perturba- 
tion theory (FT) series. The initial perturbative series 
that is obtained after the renormalization procedure is 
not the final product of the theory. This series can be 
modified and its properties can be improved on the basis 
of additional information coming from general properties 
of the quantity under consideration. In this note we con- 
sider various descriptions of the perturbative component 
in the context of an analysis of the inclusive decay of 
the r lepton. We will discuss merits and drawbacks of 
the series expansion for the i?,--ratio in terms of powers 
of the parameter as{M^) the prescription of Ref. ||] 
which uses a contour representation, and the approxima- 
tion based on the analytic approach proposed in Ref. . 

The main object in a description of the hadronic decay 
of the T-lepton and of many other physical processes is 
the correlator Il{q^) or the corresponding Adler function 
D{—q^) = —q^dIV{q^)/dq'^. The analytic properties of 
the /^-function are contained within the relation 

D{Q') = Q' - — Tz^m, (1) 

Ja (s + Q2)- 

where i?(s) = Imn(s)/7r. According to this equation, 
the Z?-function is an analytic function in the complex 
(5^-plane with a cut along the negative real axis. 

After renormalization, the perturbative expansion of 
the _D-function has the form of a power series in the ex- 



pansion parameter = 
the series has the form 



as{fJ. In the massless case 
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As is well known, this expression is unsatisfactory both 
from theoretical and practical viewpoints. Any partial 
sum of it is not renormalization-group invariant and log- 
arithms in the coefficients lead to an ill-defined behavior 
in both infrared and ultraviolet regions. 

The modification of the initial representation (|^) based 
on renormalization-group invariance reads 



D(Q2) oc 1 + d{Q^)di^o + a'(Q')d2,o 
+ d^{Q^)d3,o + -- - , 



(3) 



where d{Q^) is the running coupling. This commonly 
used modification removes some of the undesirable fea- 
tures of the expansion (Q) . A partial sum of the series (||) 
is now /^-independent. The log-terms in the coefficients of 
Eq, (||) have been summed into the running coupling and 
the series (||) can now be used in the ultraviolet region. 
However, the correct analytic properties of the partial 
sum of Eq. the principal merit of this expansion, 
are no longer valid due to unphysical singularities of the 
perturbative running coupling. 

An analytic approach proposed in Ref. |^ gives a pos- 
sible resolution of this problem. The series (||) has been 
constructed from that in Eq. (j2|) by using additional in- 
formation of a general type — the renormalization-group 
invariance of the quantity under consideration. The an- 
alytic approach makes the next logical step in the mod- 
ification of the perturbative expansion by bringing into 
consideration additional general principles of the theory 
which are reflected in the Q^-analyticity. In the frame- 
work of analytic perturbation theory (APT), the series 
has the form of a nonpower expansion 0] 



D{Q^) oc 1 + A(Q') di,o + ^2(Q') ^2,0 
+ AsiQ^) ds^o + ■ ■ ■ , 



(4) 



where Ak{Q'^) are analytic functions in the complex Q^- 
plane with a cut along the negative real axis. The expan- 
sion (0) maintains not only correct analytic properties of 



1 



the partial sum, but leads to some new remarkable fea- 
tures. For example, within the APT the renormalization 
scheme (RS) dependence of results obtained, caused by 
the inevitable truncation of the series, is reduced dras- 
tically (see details and applications to various processes 
m Rcfs. (|-[l0|). Mor cover, the analytic approach allows 
one to give a self-consistent definition of the perturbative 
expansion in the timelike region |pl]| . The APT represen- 
tation of the i?-ratio in the process of e^e~ annihilation 
into hadrons, defined in the timelike region, like the D- 
function, defined in the spacelike region, has the form of 
a nonpower expansion 



R{s) oc 1 + 2li(s) di,o + 2t2(s) d2,o 
+ 2l3(s)d3^o + -- - • 



(5) 



Both expansions, Eqs. and (||), are related term by 
term by Eq. (|l|). The functions Ak and 2lfc are close to 
each other, although they do not coincide [pTj-p^. 

The functions D{Q^) and R{s) can be written in terms 
of an effective spectral function p(a) 
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The function p{a) is calculated by taking the discontinu- 
ity of D(Q'^) across the cut. 

The leading-order expansion functions in the spacelike 
and timelike regions are 



Aim = ^ 



2li(s) = -r 

Po 



A^ 



ln(Q2/A^) A2-g2 



1 1 ln(s/A2) 
arctan 

2 TT TT 



(8) 



(9) 



The expression (g) contains the standard logarithmic 
term which coincides with the perturbative expression 
having a ghost pole at = A^. The second term in 
Eq. (||) (which appears automatically from the Kallen- 
Lehmann representation) has a non-logarithmic power 
structure, cancels the ghost pole, and provides expression 
(^ with the correct analytic properties. This term has 
a short distance origin unrelated to the operator prod- 
uct expansion and, therefore, is not inconsistent with it 
(see, for instance, the discussion in Ref. jl^). Rewritten 
in terms of the PT running coupling the second term in 
Eq. (H) has the structure of the type exp(— 1/a) and it 
therefore makes no contribution to the perturbative se- 
ries. Thus, the true analytic properties of the function 
(^ are restored by terms which are invisible in the per- 
turbation expansion, which terms, nevertheless, are felt 
by the dispersion relation, even for the spectral function 
calculated perturbatively. The regularity of the function 



(^, defined in the timelike region, has another origin. Its 
perturbative expansion contains only logarithmic terms 
proportional to powers of tt^ / In^ {s / A'^) . The Tr^-terms 
are usually incorporated into the coefhcients in the per- 
turbative expansion of R{s) without changing the expan- 
sion parameter, i.e., the form of the running coupling in 
the timelike region is taken to be the same as in the space- 
like region. In any order of that expansion the relation 
(|l|) will be violated. 

The experimentally measured i?T--ratio of hadronic to 
leptonic widths is given by 

This expression can be represented in the form of a con- 
tour integral in terms of the Adler _D-function ||l| 

(11) 

The integration in Eq. ( |l0| ) contains an interval over 
timelike momentum which extends down to small s and, 
therefore, cannot directly be calculated if one applies the 
standard parameterizations of R{s) in terms of the sin- 
gular PT running coupling. At first glance the contour 
representation ( pT| ) solves this problem because the in- 
tegration over the circle is well-defined even if the D- 
function is written in terms of as{Q^). However, in or- 
der to perform this transformation self-consistently, it is 
necessary to maintain required analytic properties of the 
ZJ-function, which are violated in the framework of stan- 
dard PT with a singular running coupling. 

Two approaches are usually used to describe the 
perturbative component in inclusive r-decay. In the 
first one called the fixed-order perturbation theory 
(FOPT), the perturbative QCD contribution 6r to Rr 
is represented in the form of a power expansion in the 
parameter a (Af^) . The three-loop approximation reads 



rFOPT 



. (M^) + Ki {Ml) -f K2 d^ (M^) , (12) 



where in the MS-scheme for three active flavors the coef- 
ficients are Ki = 5.2023 and K2 = 26.366. 

Another PT form is obtained if the perturbative ex- 
pansion (HI) is substituted into the contour integral ([Tl]). 
This gives the so-called contour-improved perturbation 
theory (CIPT) representation ||^ 



rCIPT 



^(1) (Af 2) + di ^(2) (^^2) ^ ^(3) ^M^^ (13) 



with 



(Ml) = -L / ^ (14) 
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where di = 1.640 and d2 = 6.371 in the MS-scheme. 

Both expressions, Eqs. (|l^) and ([l3|), are widely used 
in the analysis of r-decay data. However, their sta- 
tus is different. The formula (|l^ ) can be obtained self- 
consistently. In Eq. ( |Tc| ) one has to use for R{s) the initial 
perturbative approximation, like Eq. (j^) for /^-function, 
with the expansion parameter a^. Then, after integra- 
tion over s, the logarithmic terms containing ln(M^//i^) 
are removed by setting fi^ = M^. The same result is 
obtained if the contour representation ( p] ) is used and 
the D-function is taken in the form (j^) which preserves 
the required analytic properties. As for the representa- 
tion ([T^), it will be consistent with Eqs. (|lO|) and ( pT| ) if 
d{z) has analytic properties of the Kallen-Lehmann type. 
The use of the standard running coupling with unphysi- 
cal singularities in Eq. ( p^ breaks this consistency. 

In the framework of APT, where the Q^-analyticity is 
maintained, the expression for dr can be obtained both 
from the initial formula (p^ ) and from the contour repre- 
sentation (|l^) . In terms of the effective spectral function 
it has the form B 



srAPT 



(15) 
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FIG. 1. Behavior of the QCD perturbative contribution 
to Rr as a function of as{M^) (for FOPT and CIPT) or of 
aan(M^) (for APT) in the MS renormalization scheme. 

In Fig. |l|, we compare St as a function of the three-loop 
running coupling in the cases of FOPT, CIPT, and APT. 
The running coupling in the Euclidean region is defined 

by 

1 flrr 

which in lowest order coincides with Eq. (||). This plot 
corresponds to calculations in the MS-scheme. The dif- 



ference between these functions is negligible at suffi- 
ciently small as(M^) and becomes substantial with larger 
values of the coupling. The shaded area reflects the value 
Sr = 0.200 ± 0.013 ||l^ which corresponds to the exper- 
imental measurement for the non-strange channel of the 
T-data 0,|l|. 

A discussion of the perturbative contribution to in- 
clusive T-decay has been given in Ref. |p^ , in which it 
has been claimed that in order to describe the exper- 
imental data for r-decay the value of aii^{M^) should 
be taken in the range of 1.5-2.0. This conclusion con- 
tradicts our result. Fig. |^ clearly demonstrates that 
in order to reproduce the experimental data such large 
values of aan(Af^) are not required. Moreover, in the 
APT approach, the value of the analytic running cou- 
pling aan(Q^) is bounded from above and cannot exceed 
the infrared limiting value aan < t^/Po — 1-4 [||. Beyond 
this, in Ref. |19[, it was noted that this impossibly large 
value aan(M^corresponds to as(M|) ~ 0.15. How- 
ever, in order to obtain the value of as at the .Z-boson 
mass scale, the region of five active quarks should be ap- 
proached by applying a special procedure of matching 
from the three-quark region | [T^,p^ . Corresponding esti- 
mates have been given in Ref. |1C|]~ The large value of the 
scale parameter A extracted in a pure APT analysis of 
the r-data indicates that nonpcrturbative effects are not 
negligible. As it has been demonstrated in Ref. |20|, the 
light ZJ-function corresponding to the non-strange vec- 
tor channel r-data can be described by using reasonable 
effective quark masses with A 



MS 



420 MeV. 



As has been emphasized in Ref. ]19[ , the merit of the 
contour representation is that it produces expressions for 
quantities in the physical region that are not expansions 
in the parameter tt/ ln( s/A^), not a small quantity in the 
intermediate energy region. In particular, one can write 
the formula for R{s) [see Eq. (11) from Ref. [|l9|] which 
contains the expression on the right-hand side of Eq. (|^) . 
Although we agree with the authors of Ref. jl^ that it is 
important to sum up the 7r^-terms in the region of inter- 
mediate energies (see Ref. iQ) and that it is preferable 
to use an expression which sums up the 7r^-contributions 
rather than the asymptotic expression cx l/ln(s/A^), we 
note the following: The expression discussed in Ref. M] 
is simply associated with the APT timelike function ^) 
which unambiguously leads to the Euclidean function 
(^. So, we conclude that the analytic approach pro- 
vides a consistent form of expressions of the type (H) in 
the timelike region, whose evident advantage is the sum- 
mation of the TT^-contributions into a regular function, 
and results in the corresponding analytic expressions for 
the D-function in the spacelike region, where unphysi- 
cal singularities are cancelled by functions of a nonloga- 
rithmic type vanishing in a perturbative expansion. In 
other words, summation of the 7r^-contributions in the s- 
channel produces power (nonlogarithmic in Q^) contribu- 
tions in the t-channel that ensure the analytic properties 
by canceling unphysical singularities in the logarithmic 
terms. 
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A significant source of theoretical uncertainty of the re- 
sults obtained arises from the RS dependence due to the 
inevitable inclusion of only a finite number of terms in 
the PT series. There are no general principles that give 
preference to a particular renormalization scheme and 
the stability of the result has to be studied at least with 
respect to a relevant class of schemes. A virtue of the 
APT approach is its higher stability over a wide range of 
RS ijl^ . (A similar observation has been made in the con- 
text of deep- inelastic scattering sum rules In con- 
trast, results obtained in the framework of FOPT have a 
strong RS dependence. Indeed, if one performs calcula- 
tions in some scheme and then by using a RS transforma- 
tion extracts the value of a^^{M^) in the commonly used 
MS-scheme, we find that this value possesses a strong un- 
certainty. This fact is shown in Fig. ^. A point on the 
(ifi, 62)-plane corresponds to a RS, where &2 = 
is the ratio of the three-loop and one-loop coefhcients of 
the /3-function. In the shaded area there are no appro- 
priate solutions of the cubic equation (|lj) . We indicate 
the popular MS scheme and the optimal schemes: ECH, 
which recently has been applied to the r-decay , based 
on an effective charge method p^L and PMS based on a 
principle of minimal sensitivity ||23[|. A straightforward 
way of extracting the value of the coupling based solely 
on the MS scheme gives (^{M"^) = 0.337 ± 0.008. An 
analysis which takes into account the RS ambiguity leads 
to a much larger error in as(M^). The reason for the en- 
hanced RS stability of the APT scheme is its perturbative 
stability — the leading order term captures the bulk of the 
QCD correction. 
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FIG. 2. Contour plot of af^iW}) extracted from the 
T-data by using different renormalization schemes in the 
FOPT approach. 

In summary, we have made a comparative analysis of 
the merits and drawbacks of three forms of perturbative 
approximations in the context of the application of QCD 
to the inclusive decay of the r-lepton. Two of them, 
the FOPT and CIPT, use, in one way or another, the 



running coupling with unphysical singularities to param- 
eterize Rr- However, the justification of these schemes 
is rather different. The FOPT has a more solid founda- 
tion than the CIPT, which is simply self-contradictory. 
However, due to the strong RS dependence, the FOPT 
leads to a large uncertainty in extracting QCD param- 
eters from the data. We have put forward arguments 
in favor of APT. Calculations in the framework of this 
approach are self-consistent and considerably reduce the 
theoretical uncertainty of the results obtained; this is im- 
portant in extracting the nonperturbative component of 
the QCD description of the process. 
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